home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / SGECO.z / SGECO
Text File  |  1996-03-14  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. SSSSGGGGEEEECCCCOOOO((((3333FFFF))))                                                            SSSSGGGGEEEECCCCOOOO((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      SGECO   - SGECO factors a real matrix by Gaussian elimination and
  10.      estimates the condition of the matrix.
  11.  
  12.      If  RCOND  is not needed, SGEFA is slightly faster.  To solve  A*X = B ,
  13.      follow SGECO by SGESL.  To compute  INVERSE(A)*C , follow SGECO by SGESL.
  14.      To compute  DETERMINANT(A) , follow SGECO by SGEDI.  To compute
  15.      INVERSE(A) , follow SGECO by SGEDI.
  16.  
  17.  
  18. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  19.       SUBROUTINE SGECO(A,LDA,N,IPVT,RCOND,Z)
  20.  
  21. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  22.      On Entry
  23.  
  24.      AAAA REAL(LDA, N)
  25.         the matrix to be factored.
  26.  
  27.      LLLLDDDDAAAA INTEGER
  28.         the leading dimension of the array  A .
  29.  
  30.      NNNN INTEGER
  31.         the order of the matrix  A .  On Return
  32.  
  33.      AAAA an upper triangular matrix and the multipliers
  34.         which were used to obtain it.
  35.         The factorization can be written  A = L*U , where
  36.         L  is a product of permutation and unit lower
  37.         triangular matrices and  U  is upper triangular.
  38.  
  39.      IIIIPPPPVVVVTTTT INTEGER(N)
  40.         an integer vector of pivot indices.
  41.  
  42.      RRRRCCCCOOOONNNNDDDD REAL
  43.         an estimate of the reciprocal condition of  A .
  44.         For the system  A*X = B , relative perturbations
  45.         in  A  and  B  of size  EPSILON  may cause
  46.         relative perturbations in  X  of size  EPSILON/RCOND .
  47.         If  RCOND  is so small that the logical expression
  48.         1.0 + RCOND .EQ. 1.0
  49.         is true, then  A  may be singular to working
  50.         precision.  In particular,  RCOND  is zero  if
  51.         exact singularity is detected or the estimate
  52.         underflows.
  53.  
  54.      ZZZZ REAL(N)
  55.         a work vector whose contents are usually unimportant.
  56.         If  A  is close to a singular matrix, then  Z  is
  57.         an approximate null vector in the sense that
  58.         NORM(A*Z) = RCOND*NORM(A)*NORM(Z) .  LINPACK.  This version dated
  59.      08/14/78 .  Cleve Moler, University of New Mexico, Argonne National Lab.
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. SSSSGGGGEEEECCCCOOOO((((3333FFFF))))                                                            SSSSGGGGEEEECCCCOOOO((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      Subroutines and Functions LINPACK SGEFA BLAS SAXPY,SDOT,SSCAL,SASUM
  75.      Fortran ABS,AMAX1,SIGN
  76.  
  77.  
  78.  
  79.  
  80.  
  81.  
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.